Strength of Correlations in electron and hole doped cuprates 
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High temperature superconductivity was 
achieved by introducing holes in a parent com- 
pound consisting of copper oxide layers separated 
by spacer layers. Realizations of this phenomena 
has been achieved in multiple crystal structures 
and has been the subject of numerous investi- 
gations and extensive controversy. In a small 
number of copper oxide based materials, it is pos- 
sible to dope the parent compound with electrons 
[TI43J , and their physical properties are bearing 
some similarities but also significant differences 
from the hole doped counterparts. For example, 
in the electron doped materials the antiferro- 
magnetic phase is much more robust than the 
superconducting phase, while the normal state 
has a resistivity with a quadratic temperature 
dependence which is expected in normal metals 
rather than the linear temperature dependence 
observed in the hole doped systems. Here, we use 
a recently developed first principles method, to 
study the electron doped cuprates and elucidate 
the deep physical reasons why their behavior is 
so different than the hole doped materials. The 
crystal structure of the electron doped materials, 
characterized by a lack of oxygen in the apical 
position, results in a parent compound which is 
a Slater insulator. Namely, a material where the 
insulating behavior is the result of the presence 
of magnetic long range order. This is in sharp 
contrast with the hole doped materials, where 
the parent compound is a Mott charge transfer 
insulator, namely a material which is insulating 
due to the strong electronic correlations but not 
the magnetic order. We study the evolution of 
the angle resolved photoemission spectra and 
the optical properties of the normal state of 
the electron doped cuprates as a function of 
doping, clarifying how their unique position 
close to, but below the metal to charge transfer 
insulator transition, accounts for their surprising 
differences from the hole doped cuprates. 

We study the prototypical electron doped compound 
Nd2-a^Cea^Cu04 (NCCO). This is a single-layer material 
and therefore we compare it to La2-ajSra;Cu04 (LSCO), 
a single layer hole doped material in the related T struc- 
ture. We use a realistic theoretical approach, e.g. the 
Local Density Approximation combined with the Dy- 
namical Mean Field Theory (LDA+DMFT) [5l|6], and 
focus on intermediate energy quantities. The results for 
the electronic spectral functions A{uj) = — l/7rImG(co') 



of NCCO are shown in Fig. [T] The DMFT equations 
have two solutions shown in panels a,b. The first one 
is non-magnetic and the second is antiferromagnetically 
ordered. Notice that the ordered solution, which is stable 
up to temperature 1400 K, has a charge transfer gap of 
1.2eV. When comparing to experiment, it is important 
to keep in mind that two dimensional compounds are 
not able to sustain infinite-range magnetic order at a fi- 
nite temperature Therefore, the Neel temperature within 
DMFT should be interpreted as the temperature below 
which the magnetic correlations become long but remain 
finite. This temperature can be much higher then the ac- 
tual Neel temperature of the material, which is controlled 
by the magnetic exchange between the two dimensional 
copper oxide layers; and vanishes for a well separated 
copper oxide planes. DMFT also allow us to investigate 
the underlying paramagnetic solution, which describes a 
material in the absence of long range magnetic order. 
This solution is metallic, and hence the magnetic long 
range order is responsible for the insulating nature of the 
compound, as surmised by Slater [8 . It is worth noting 
that the oxygen orbitals carry no magnetic moment be- 
cause of the polarization with the two copper neighbors, 
which have the opposite moments. 

In Fig. ^ we resolve the low energy part of panel b 
in momentum space. We find two dispersive coherent 
peaks separated by the charge transfer gap. In Fig. [l]i 
we resolve the low energy part of panel b in orbital space. 
We find that the occupied and empty states have mostly 
copper d character with significant oxygen p admixture. 
The top of the lower band occurs at (7r/2, 7r/2), while the 
bottom of upper band apperars M = (tt, 0), therefore the 
gap is indirect (see yellow arrow in panel c). Those two 
bands can also be obtained in the simpler Hartree Fock 
approximation, though the size of the gap is overesti- 
mated in a static mean-field. 

In Figs, [l^ and f we show the electronic structure at 
10% electron doping. The compound is still magnetic. 
Some aspects of the doped electronic structure can be 
understood in terms of the Hartree Fock rigid band pic- 
ture, the holes appear first in the M = (tt, 0) point, there 
are additional narrow quasiparticle states close to but 
below the Fermi level. The two peak structure of the 
low energy quasiparticles is clearly seen in the orbital 
resolved angle integrated spectra of Fig{T] f. 

Notice also the presence of the pseudogap around 
(7r/2, 7r/2), that is a signature of the presence of magnetic 
long range order, and is also observed in experiments [7 
(Fig. [l^, right panel), and compares well to our theoret- 
ical calculations (Fig. [It, left panel). The pseudo-gap 
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FIG. 1: Density of electronic states in NCCO : (Col- 
ors online) Spectral function of NCCO a) nonmagnetic and b) 
magnetic) at zero doping in the extended energy range. There 
is a gap of 1.2eV when magnetism is allowed. We clearly ob- 
serve the presence of the lower (upper) Hubbard bands LHB 
(UHB), and of a quasi-particle peak in the non-magnetic com- 
pound (QP). Frequency dependant spectral weight A(k,w) 
along F-M-K-F obtained by LDA+DMFT at c) integer- 
filling and e) 10% electron doping (NCCO). The partial den- 
sity of states of the d and p orbitals are shown on the right 
panels d) and f). The lines in c) and e) are the LDA rigid 
bands, plotted in the folded Brillouin Zone, for comparison 
with DMFT. The corresponding integrated weight is plotted 
in d) and f). g) side by side comparison of A(k,w), along 
the path F — i^ as depicted in the inset, obtained theoreti- 
cally (left side) and experimentally from Ref. [7] (right side) 
at 13% doping, h) comparison of A(w) for a fixed k point 
k — (37r/4, 7r/4) (shown in the inset) obtained theoretically 
(lower curve) and experimentally from Ref. 7 (upper curve). 
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FIG. 2: Optical properties of NCCO : (Colors online) 



a) Theoretical optical conductivity in Q 



for NCCO 



doped with charge carries concentration at integer filling (red 
curve) and at 10 percent doping (blue curve). At integer- 
filling, the system is gapped and we observe an optical gap of 
about 1.5eV, that is larger than the charge gap ~ 1.2eV. The 
left inset is a magnifying glass on the data. For comparison we 
also show (right inset) the optical conductivity obtained by 
infrared optics ^9 . b) We show the dimensionless integrated 
optical conductivity Neff for LDA+DMFT done on LSCO 
and NCCO, obtained with a cutoff A = 1.2 (A = 1.5) for 
NCCO (LSCO). Experimental data for LSCO [TO] (red circles) 
and NCCO [11 (open diamonds) are shown. 



closes in the paramagnetic compound. In Fig. [T)i, we 
also show the spectral weight resolved in frequency for a 
fixed k = (37r/4, 7r/4) point (lower curve). We emphasize 
the ability of DMFT to resolve the multiple peak struc- 
ture recently measured in experiments [7 (Fig. [iji upper 
curve), that cannot be described in a mean-field picture 
like Hartree Fock. 

We now turn to the theoretical optical conductivity, 
displayed in Fig|2^, which is in good qualitative agree- 
ment with the experimental results reproduced in the 
right inset of Fig [2^. The undoped compound has a 
sharp onset at an energy of the order of 1.5 eV which we 



interpret as the direct gap (slightly larger than the charge 
transfer gap in Fig. [T]i and a tail resulting from indirect 
transtions. Doping introduces several new features. The 
1.5eV optical peak loses weight which is transferred to 
lower energy in the form of a Drude peak and a mid in- 
frared peak at cj ~ {).2eV . 

This peak can be explained in our theory by the pres- 
ence of magnetism, that leads to a reconstruction of the 
quasi-particle band structure as shown in Fig. fib), in 
agreement with Ref. [12 , and additionally our study al- 
low to connect this peak with the spectral weight below 
the Fermi energy at the M = (tt, 0) point. 

Finally, the optical conductivity display a peak in the 
magnetic solution at a much lower frequency uj ^ 0.03eV 
(see left inset of Fig. [2^), which is connected to the 
pseudo-gap at k = (7r/2,7r/2), and is also observed in 
experiments [11 . 

To quantify the rate of the redistribution of op- 
tical spectral weight, is useful to consider the effec- 
tive electron number per Cu atom defined by N^rr = 

^^^ /q cr' {uj)duj^ where rUe is the free electorn mass, 
and V is the cell Volume containing one formula unit. 
Neff is proportional to the number of electrons involved 
in the optical excitations up to the cutoff A. Our results 
for (electron doped) NCCO are displayed in the left hand 
side of Fig. ^p and are compared to experimental data 
taken from Ref. pjj. Notice a favorable agreement be- 
tween the theory and experiment, for which the use of the 
realistic electronic structure is essential. In FigJ2J3 (right) 
we also show our theoretical optical spectral weight for 
a prototypical hole doped cuprate (LSCO), computed in 
Ref. [4], which is a Mott type insulator, and we compare 
it to experimental data [TO^ . The results for a Slater type 
insulator NCCO and Mott type insulator LSCO are very 
similar, hence we can conclude that the optical spectral 
weight is not a precise measure of the strength of corre- 
lations. 

In a Slater picture, the onset of antiferromagnetism 
reduces the Coulomb correlations (double occupancy) at 
the expense of the kinetic energy. The opposite is true in 
a Mott insulator. Consequently, in a Slater insulator the 
kinetic energy becomes less negative as the temperature 
decreases while the opposite happens in a Mott insulator. 
The kinetic energy as a function of temperature is readily 
available in the theory and is displayed in Fig [sk (top 
panel). 

A closely related quantity A^e// is experimentally ac- 
cessible and has been measured for this compound. In 
Fig [3J3 (top panel) we plot both the experimental [TT] 
and our theretical data, which are in very good agree- 
ment. Experiments of Ref. [TT] confirm that A'e// is 
closely related to minus the kinetic energy. It has there- 
fore already been confirmed that NCCO is below the 
metal to charge transfer insulator transition in experi- 
ments. 
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FIG. 3: Kinetic energy and magnetic phase dia- 
gram: (Colors online) a) Theoretical kinetic energy vari- 
ation —AEkin = — E kin (T) -\- E kin {S9° K) of the electron (top 
panel) and the hole (bottom panel) doped compounds with 
respect to temperature, b) Normalized variation of A'e//, 
ANeff = Neff{T) - Ne ff {89'' K) of the electron (top panel) 
and the hole (bottom panel) doped compounds with respect 
to temperature. The dashed line are extracted from experi- 
ments done on NCCO [TT]. c) Mean- value of the staggered 
magnetization obtained for LSCO and NCCO by both single 
site (DMFT) and cluster DMFT (c-DMFT). Experimental 
values M(6)/Mo (Mo = M(S = 0)) for NCCO 13^ and LSCO 
[14] are also shown, and for comparison with DMFT, we as- 
sume Mo = Mdmft{S = 0), where Mdmft is the magnetic 
moment at doping obtained by single site DMFT. We also 
show the magnetization obtained in the static hartree-fock 
approximation (dashed line) are shown for comparison. 



In Fig. [3^ (bottom panel) we also plot the tempera- 
ture dependence of the kinetic energy and A^e// for the 
hole doped Mott insulator (LSCO). We find the oppo- 
site trend compared to the lightly electron doped NCCO 
compound, confirming that LSCO is a doped Mott insu- 
lator. 

Finally, in order to check the validity of our single site 
DMFT approach, we computed the magnetic phase dia- 
gram and magnetic moment within the single site and 



two-site cluster DMFT (see Fig. [3Jc). Our data for 
NCCO are displayed on the left side of Fig. |3]c) and are 
compared to experimental data from Ref . [13] . The agree- 
ment is very good for both approaches, hence justifying 
the use of the single sited DMFT approach. 

The right site of Fig. [3Jc) displays theoretical and ex- 
perimental magnetic moment [14 for hole doped LSCO. 
We found significant difference in the region of stability 
of the magnetic state between the single site and clus- 
ter DMFT {S < 10%) , hence the dynamical short range 
correlations - absent in single site approach - are very 
important for LSCO in the underdoped regime, but not 
for NCCO. 

We also carried out the Hartree-Fock calculation, and 
we found that in this static approach the magnetization 
vanishes only at unrealistic large doping 6 ^ 50% for both 
NCCO and LSCO (dashed line in Fig.jsjc), which points 
towards the important role of dynamic correlations at 
finite electron and hole doping. 

Our theory sheds light on many puzzling observa- 
tions. It was noticed within the context of the one 
band model that the doping dependent Hubbard inter- 
action was needed to reproduce experiments [15 . In our 
view, this is a description, within a one band model, of 
the rapid metallization process in a Slater-like insula- 
tor. A realistic LDA+DMFT treatment of the multiband 
Hamiltonian does not require doping dependent interac- 
tions. The strong sensitivity of the electronic proper- 
ties to magnetism, which results from the fact that the 
physics of NCCO is closer to that of a Slater insulator 
than of a Mott insulator, gives rise to additional experi- 
mental predictions. We expect that at high temperatures 
(of the order of the DFMT mean field temperatures) 
there should be appreciable transfer of spectral weight 
from high energies to low energies with decreasing tem- 
perature, as the magnetic correlations are weakend, and 
a substantial decrease of the Mott Hubbard gap and the 
weights in the Hubbard bands. This can be checked by 
extending transport and photoemission studies in lighlty 
doped electron doped cuprates from 600 to 1100 K. 

In materials where the strength of the interaction is 
below the critical value need to produce a Mott insulat- 
ing state, the effective interactions can renormalize down, 
allowing superconductivity and magnetism to coexist mi- 
croscopically rather than exclude each other. This coex- 
istance provides a good description of the Raman scatter- 
ing of the electron doped superconductors [19 . Finally, 
in some NCCO films the parent compound was found to 
be metallic. It would be interesting to study the mag- 
netic properties of these materials to see if this metallic- 
ity correlates with a substantial decrease of the magnetic 
correlations [20 . 
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METHOD 

The LDA calculation was carried out with the PWSCF 
package [21 , which employs a plane- wave basis set and 
ultrasoft pseudopotentials ^. We downfold to a three 
band model, containing copper dx2-y2 and two oxygen p^r 
orbitals using the maximally localized Wannier functions 
method (MLWF) [23^. The results of this procedure can 
be summarized in the following three band Hamiltonian: 



n 



J2 ^r/4aC.-/3a + epNp^ 

ija,{a,l3)e{Px,Py,dx2-y2) 

{ea - E"-) Nd + UdJ2 f^idtnidi (1) 



where i and j label the Cu02 unit cells of the lattice, 
and t^j are the hopping matrix elements, and e^^, Cp are 
the on-site energies of the d and p orbitals. The onsite 
Coulomb repulsion U on the dx2-y2 orbital was fixed to 
Ud = SeV. Finally, we note that e^ — €p plays the role of 
an effective onsite repulsion /7 in a Hubbard model pic- 
ture. The LDA+DMFT method, accounts for the cor- 
relations which are included in both LDA and DMFT 
by a double counting correction to the d-orbital energy, 
Edc which amounts to a shift of the relative positions 
of the d and p levels which we take to be doping inde- 
pendent and Edc = 4.8eV (3.12eV) for NCCO (LSCO). 
The LDA downfolded hopping parameters are presented 
in the online material. The Green function of the three 
band model is given by: 



G]^{iuJn) = {i^n + /i - Hk - T^{iuJn)y 



(2) 



where Hk is the Fourier transform of the 1-Lt in Eq. ([T]) 
and X) is the DMFT self-energy matrix, which is assumed 
to be local, and nonzero only for the d orbital. The self 
energy in Eq. (pi is obtained by solving an Anderson im- 
purity model subject to the DMFT self-consistency con- 
dition 



■ Ej, 



^{icu) - A(zcj) Nk 



^ ^ Gf (zo;), (3) 



kGBZ 



where the sum runs on the first Brillouin Zone (BZ). The 
self-energy I] is obtained by solving an Anderson im- 
purity using the continuous time quantum Monte Carlo 
impurity solver [24l [25] (CTQMC) and the exact diago- 
nalization solver [5] (ED). Calculations have been carried 



out at temperature T = 89° K, when not specified in the 
text. Finally, the optical conductivity is given by 



^ ^ Nk^ he J 



f{x-uj) - f{x) 



X Tr (pka {x - cj) VkPka (a^) Vk) , (4) 

where c is the interlayer distance, the density matrix p is 
P\ca{x) = 2^7 (C^L(^) ~ Gka(^)j, and the bare vertex 
is obtained by following the steps of Ref ^6^ . 
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